Dynamical Localization of Bose-Einstein condensate in Optomechanics 
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We explain dynamical localization of Bose-Einstein condensate (BEC) in optomechanics both in 
position and in momentum space. The experimentally realizable optomechanical system is a Fabry- 
Perot cavity with one moving end mirror and driven by single mode standing field. In our study we 
analyze variations in modulation strength and effective Plank's constants. Keeping in view present 
day experimental advances we provide set of parameters to observe the phenomenon in laboratory. 
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I. INTRODUCTION 



Cavity-optomechanics deals with the interaction of an 
optical field in a resonator with confining mirrors [1]. Re- 
cent experimental advances make it possible to couple 
cold atoms and Bose-Einstein condensates (BEC), me- 
chanical membrane and nano-sphere with the optome- 
chanical systems [2]. Hence, these hybrid-opto mechan- 
ical systems are playground to study phenomena re- 
lated to mirror-field interaction, and atom-field interac- 
tion which provide founding principles to develop nu- 
merous sensors, and devices in quantum metrology. In 
opto-mechanics the mechanical effects of light lead to 
cool the motion of a movable mirror to its quantum me- 
chanical ground state [3-5], and to observe strong cou- 
pling effects [6-8]. In earlier work optomechanical sys- 
tems were suggested to develop gravitational wave de- 
tectors [9] and to measure displacement with large accu- 
racy [10], and in recent research work these systems lead 
to develop optomechanical crystals [11]. Recent theo- 
retical discussions and simulations on bistable behavior 
of BEC-optomechanical system [12], high fidelity state 
transfer [13, 14], steady-state entanglement of BEC and 
moving end mirror [15], macroscopic tunneling of an op- 
tomechanical membrane [16] and role reversal between 
matter- wave and quantized light field are guiding towards 
new avenues in cavity optomechanics. In this paper, we 
discuss the dynamical localization of ultra-cold atoms or 
dilute BEC trapped inside a single mode driven optome- 
chanical cavity with one fixed mirror and the other as 
a movable mirror. Dynamical localization is an inter- 
esting phenomenon in periodically driven nonlinear sys- 
tems [18-22] and provides quantum mechanical limits on 
classical diffusion of wave packet. We show that for dif- 
ferent modulation regimes the dynamical localization is 
observed both in position and in momentum space. The 
phenomenon of dynamical localization emerges as the hy- 
brid optomechanical system is explicitly time dependent, 
hence, the single mode laser field provides spatially peri- 
odic potential to BEC with phase modulation[17] due to 



modulated end mirror. 

In the paper, the model of the system is presented in 
section II. In section III we derive the Langevin equa- 
tion and in section IV, we explain dynamical localization 
of the condensates in the system. In section V, we ex- 
plain dynamical localization as a function of modulation 
amplitudes. Results are summarized in sect ion- VI. 



II. THE MODEL 

We consider a Fabry-Perot cavity of length L with a 
moving end mirror driven by a single mode optical field 
of frequency uo p and BEC with N-two level atoms trapped 
in an optical lattice potential [23, 24]. Moving end mirror 
has harmonic oscillations with frequency uo m and exhibits 
Brownian motion in the absence of coupling with radia- 
tion pressure. 

The Hamiltonian of BEC-optomechanical system is, 



H = Hrr 



(1) 
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where, H m describes the intra-cavity field and its cou- 
pling to the moving end mirror, H a accounts for the 
BEC and its coupling with intra-cavity field while, Ht 
accounts for noises and damping associated with the sys- 
tem. The Hamiltonian H m is given as [25], 

Hm = HA c c^c + ^ (p 2 + <f ) - ^cq - itvq{c - c f ) , (2) 



where, first term is free energy of the field, A c = uo c — uo p 
is detuning, here, uj c is cavity frequency and (c) are 
creation (annihilation) operators for intra-cavity field in- 
teracting with mirror and condensates and their com- 
mutation relation is [c, c^] = 1. Second term represents 
energy of moving end mirror. Here q and p are di- 
mensionless position and momentum operators for mov- 
ing end mirror, such that, [q,p] = i. Intra-cavity field 
couple BEC and moving end mirror through radiation 
pressure. Third term represents such coupled energy of 
moving end mirror with field and £ = \/2(uj c / L)x$ is 
the coupling strength where, xq — ^h/2mcj m , is zero 
point motion of mechanical mirror having mass m. Last 
term gives relation of intra-cavity field and output power 



2 



1 77 1 = ^/ Pn/fajOp, where, P is the input laser power and 
k is cavity decay rate associated with outgoing modes. 

Now the Hamiltonian for BEC and intra-cavity field 
and their coupling is derived by considering quantized 
motion of atoms along the cavity axis in one dimensional 
model. We assume that BEC is dilute enough and many 
body interaction effects are ignored. We have 



fid 2 



2m a dx 2 



HUqc^ ccos 2 kx I ip(x)dx, 



(3) 

here, ^(V^) is annihilation (creation) operator for bosonic 
particles and Uo = 9o/A a is the far off-resonant vacuum 
Rabi frequency. Here, A a is far-off detuning between 
field frequency and c^o, atomic transition frequency, m a is 
mass of an atom, and k is the wave number. Due to field 
interaction with BEC, photon recoil takes place that gen- 
erates symmetric momentum ±2lhk side modes, where, I 
is an integer. We assume that field is weak enough which 
causes low photon coupling, therefore only th and I s * 
order modes are excited while, higher order modes are 
ignored. Now i[j is depending upon these two modes [12] 
and defined as, ip(x) = d J r \/2cos{2kx)b/ L. Where, a and 
b are annihilation operators for th and ±2Hk th modes re- 
spectively and are related as a)d + tfb = TV, where, N 
is the number of bosonic particles. As population in th 
mode is much larger than the population in 1 st order 
side mode, we write a) a ~ N or a and a) — >• y/~N . This is 
possible when side modes are weak enough and for that 
matter can be ignored. 

By using ^(x) defined above in Hamiltonian, we write 
the Hamiltonian governing the field-condensate interac- 
tion as, 



^etc+^(P 2 + g 2 )+^ctcQ (4) 



here, first term describes energy of field due to the con- 
densate. We also assume large atom-field detuning A a , 
therefore, excited atomic levels are adiabatically elim- 
inated. Second term expresses the energy of the con- 
densate in the cavity following harmonic motion. Here, 
P = ^75 — and Q = ~^(b + ¥) are dimension- 
less momentum and position operators for condensate 
which are related as [Q,P] = i and Q = 4u r where, 
uj r = hk 2 /2m a is recoil frequency of an atom. Last term 
in eq. (4) describes coupled energy of field and con- 
densate with coupling strength £ sm = ^^/h/m srn Auj r , 
where, m sm = Hu 2 /(L 2 NUQUJ r ) is the effective mass of 
side mode. 



III. LANGEVIN EQUATIONS 

The Hamiltonian Ht accounts for the effects of dissi- 
pation in the intra-cavity field, damping of moving end 
mirror and depletion of BEC in the system via standard 
quantum noise operators [30]. The total Hamiltonian 



H leads to coupled quantum Langevin equations for po- 
sition and momentum of moving end mirror and BEC, 
viz., 

dc 

dt 
dp 

~dt 
dq 

~dt 

dP 

~dt 



(iA + i^q - i^sm - k)c + V + V2ko~, 

-Umq - - JmP + fB, 



4:UJ r Q (sm^ C ^ srn P -\- firm 



£-4 



4uj r P 
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2m- 



(5) 



In above equations A = A c — NUq/2, whereas a\ n is 
Markovian input noise of the cavity field. The term j m 
gives mechanical energy decay rate of the moving end 
mirror and /# is Brownian noise operator [31]. The term 
7sm represents damping of BEC due to harmonic trap- 
ping potential which effects momentum side modes while, 
fiM and fiM are the associated noise operators assumed 
to be Markovian. 

The coupled equations of motion which govern the dy- 
namics of the moving end-mirror and the BEC are quan- 
tum Langevin equations (5). In adiabatic approximation 
we consider no thermal excitation, for the reason, we ig- 
nore noises and radiation effects resulting from optical 
damping [1, 24, 32]. Therefore, in short time limit we 
ignore the mechanical damping, as a result the coupled 
equations of motion for the moving end-mirror and BEC 
reveal their coupled nonlinear oscillator dynamics, that 
is, 



d 2 q 2 . 



d 2 Q 
dt 2 



-4uj 2 Q 



4,UJr^smV 



(A + ^-^mQ) 2 



(6) 



The corresponding effective Hamiltonian is described as 
H e ff = K + V, where, K = p 2 /2, here, k = 1 is dimen- 
sionless Planck's constant. The potential V is governed 
by Eq. (7). We assume that for a weak coupling the 
moving end mirror behaves like harmonic oscillator with 
frequency (jj m and q = qo cos(cj m t), where, qo is the max- 
imum displacement from mean position. We introduce 
some dimensionless parameters defined as, 7 = 4o; r /o; m , 
/3 = if/n 2 , 11 = A/k, iii = ism/K, A = ^f-qo and di- 
mensionless time, 4r = oj m t. Hence, from Eq. (7) we get, 



dr 2 



l + [^-/zi{Q-Acos(4r)}] 2 



dQ. 

(8) 



Now, we write the Hamiltonian as, 

*■ 1 r 47&3/U, 



[M-Mi{Q-Acos(4r)}] 2 
(9) 



dQ. 
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By using transformation X = Q — Acos(4r), we find an 
effective Hamiltonian for the condensate, that is, 

H e ff — -P 2 -\--X 2 +X\ e ff cos(4r) — arctan(/i— /iiX), 
2 2 fi\ 

_ (10) 

where, 7 m = and X e ff = (1 + ^|)A is the effective 
modulation. 

In our later work we consider the power of external field 
P = 0.0164raW with frequency uj p = 3.8 x 2tt x 10 14 Hz 
and wave length X p = 780nra. Coupling of external field 
and intra-cavity field is r] = 18.4 x 2ttMHz and frequency 
of intra-cavity field is considered uo c = 15.3 x 2tt x 10 14 i7z 
which produces recoil of uo r = 3.8 x 2itkHz in atoms 
placed in cavity of length I = 1.25 x 10 _4 m and having de- 
cay rate ft = 1.3 x 2ttMHz. Moving end mirror of cavity 
should be perfect reflector that oscillates with a frequency 
cu m = 15.2 x 2irkHz. The mirror-field and condensate- 
field coupling strengths are, respectively, £ = 14.39KHz 
and £ sm = 15.07MHz. Detuning of the system is taken 
as A = A c + ^f- 0.52 x 2ttMHz, where vacuum 
Rabi frequency of the system is Uo = 3.1 x 2ttMHz 
and number of ultra cold atoms placed in the BEC- 
optomechanical system are N = 2.8 x 10 4 [2, 26-29]. 



IV. DYNAMICAL LOCALIZATION OF THE 
CONDENSATE 

The quantum dynamics of dilute BEC is explored by 
evolving a Gaussian material wave packet, ^(X), defined 
at r = with, AX as its initial dispersion in the position 
space and AP corresponding to its initial dispersion in 
the momentum space. 

Dynamics of the condensate is shown in Fig-1. Here, 
Fig- 1(a) shows classical and quantum dispersion in mo- 
mentum space as a function of time. Initially, the clas- 
sical and quantum dispersion in momentum space fol- 
lows the classical diffusion law but after quantum break 
time, quantum dispersion of the material wave packet 
saturates and oscillates around a mean value while, clas- 
sical dispersion increases continuously with time, follow- 
ing t a law, where a = 0.61 showing anomalous diffusion 
in classical domain. Similarly, from Fig- 1(b), we note 
that quantum dispersion in position space saturates af- 
ter quantum break time and fluctuates around average 
value while, classical dispersion in position space that 
shows anomalous diffusion as well. 

Here, Fig.-l(d,e) show classical and quantum time av- 
eraged probability distributions of BEC both in momen- 
tum and in position coordinates, respectively, for mod- 
ulation amplitude A e // = 1.05 x 10 _5 m. We calculate 
probabilities in momentum coordinate W(P) = |^(P)| 2 
and position coordinate W(X) = |\I/(X)| 2 after an evolu- 
tion time t = 4.19 x 10 -2 sec. From the numerical results, 
we observe that the classical probability distribution in 
momentum and in position space are broader than their 
quantum mechanical counterpart. In contrast quantum 



mechanical probability distributions are maximum at the 
points where the matter wave was initially placed and as 
we traverse the regions away from the initial mean val- 
ues the probability of finding the matter wave both in 
momentum and in position space decreases. The mixed 
phase space with chaotic and regular regions appear due 
to the presence of modulated end-mirror that modifies 
the evolution of the condensate in effective nonlinear po- 
tential seen by the condensate in the cavity, as shown in 
Fig-l(f). 

As probability distributions are marginal distributions 
of phase space, in quantum mechanical probability distri- 
butions the peak positions both in momentum and in po- 
sition coordinates correspond to the underlying nonlinear 
resonances in the Poincare surface of section. Maximum 
fraction of the condensate is trapped with momentum in 
the range —3 x 10 _2T to 3 x 10~ 27 kg mj 'sec and in the 
position space —3 x 10 _5 m to 3 x 10 _5 m because of reg- 
ular regions in this domain as shown in Fig- 1(f). In mo- 
mentum distribution, peak at P ~ 1.4 x 10 -30 /cg m/ sec 
corresponds to the regular region at same momentum 
value. Similarly, in position space distribution peaks at 
X ~ -6.2 x 10" 6 and X ~ 6.5 x 10" 6 corresponds to the 
regular regions appearing in phase space. 

The part of the wave packet which is not in regular re- 
gion experiences exponential decay which indicates that 
the probability of finding the condensate decreases ex- 
ponentially in the regions away from the classical reso- 
nances. Maximum probability of finding the condensate 
is near the primary resonance as it originates initially 
from there. The exponential decay of distribution is one 
of the fundamental evidence of dynamical localization. 

Fig-1 (c) shows mean quantum mechanical energy < 
H > of the condensate as a function of time for A e // = 
1.05 x 10 _5 m. The mean energy of the condensate shows 
very small fluctuation around a constant value and shows 
quantum recurrences as a function of time. 

The consequence of dynamical localization in position 
and momentum space becomes prominent as we study 
the spatio-temporal dynamics of the condensate, pre- 
sented in Fig- 2. We show space-time dynamics in mo- 
mentum space both in quantum mechanical and classical 
space, respectively, in Fig. 2(a) and 2(b). Whereas, we 
plot space-time dynamics in position space both quantum 
mechanically and classically, respectively, in Fig. 2(c) and 
2(d). As the time evolves classical distributions both 
in momentum space and in position space melt down 
and spread over the entire available space with maximum 
probability in the resonances. On the other hand, quan- 
tum mechanical probability distributions in momentum 
space and in position space initially follow classical be- 
havior, however, beyond quantum break time localization 
limits the spread. The maximum probability of finding 
the material wave packet stays around the nonlinear res- 
onances. The material wave packet tunnels to other res- 
onances which provide the exponential decay as shown 
in Fig. l(d,e), therefore probability of finding the con- 
densate in neighboring resonances increases with time. 
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FIG. 1: (a) shows dispersion Ap in momentum space and (b) shows dispersion Ax in position space vs time for fixed modulation 
and gives comparison in quantum mechanical and classical behavior, (c) Mean energy of the condensate as a function of time, 
(d) Classical and quantum mechanical probability distribution in momentum space, (e) Classical and quantum mechanical 
probability distribution in position space of condensate for fixed time t = 4.19 x 10 -2 sec. (f) Phase Space have initial 
distribution from -4 to 4 both in position space P x 10 -26 kgm/ sec and momentum space X x 10 _5 m. In all these figures 



modulation A 



eff 



1.05 x 10" 5 m. Moreover, \i = -0.4, fn = 2, ft = 1.8 and 7^ 



0.6034. 



V. EFFECTS OF MODULATION 




-3 „ 3-3 3 



P (xlO^Kgrn/sec) X(x10 5 m) 

FIG. 2: Classical and quantum mechanical probability distri- 
bution versus time, (a) Quantum mechanical distribution in 
momentum space and (b) in position space, (c) Classical dis- 
tribution in momentum space and (d) in position space. The 
modulation is fixed as A e // = 1.05 x 10 _5 m. The remaining 
parameters are same as in Fig-1. 



We numerically study the classical dispersion and 
quantum dispersion of the condensate both in momen- 
tum space Fig-3(a) and in position space Fig-3(b) as a 
function of modulation strength for a fixed time. For 
a small value of modulation classical and quantum dis- 
persions deviate slightly both in position and in momen- 
tum space. As we increase the modulation in the sys- 
tem, momentum and position dispersions show contrast- 
ing behavior in classical and in quantum mechanical do- 
main. In classical dynamics dispersion both in position 
space and in momentum space increases while, quantum 
dispersion displays saturation. The quantum dispersion 
manifests periodic oscillations in momentum space and 
in position space. At the minima of these oscillations 
we observe strong localization at particular modulation 
values where the difference between classical dispersion 
and quantum mechanical dispersion is maximum. As 
the value of rescaled Planck's constant, % reduces, the 
time evolution of the quantum dispersion both in posi- 
tion and momentum space approaches to classical disper- 
sion and periodic oscillation in the quantum mechanical 
dispersion vanish as well, hence, the periodic behavior is 
strictly quantum phenomena which may display quantum 
revivals in dynamical systems. 
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(x10*m) 2.6 



FIG. 3: Left figure shows classical and quantum dispersion in 
momentum space and right figure shows classical and quan- 
tum dispersion in position space as function of modulation. 
These dispersion in momentum and position space are ob- 



served after an evolution time t — 
parameters are same as in Fig. 1. 



2.09 x 10~ 2 sec. The other 



coupled to the moving end mirror through the cavity field 
and displays dynamically localization both in position 
and momentum space. A contrast between classical and 
quantum dynamics is the evidence of the existence of dy- 
namical localization and its presence is the signature of 
quantum chaos in a dynamical system. This behavior 
is different from the dynamical localization of ultra-cold 
atoms in modulated optical field where localization is ob- 
served only in momentum space [20]. Spatio-temporal 
dynamics of condensate in position and momentum space 
confirm our finding as well. We have also calculated the 
dispersion in position and momentum space as a function 
of modulation for fixed time. The dynamical localization 
phenomenon is realizable experimentally in position and 
in momentum space by using presently available experi- 
mental parameters. 



VI. CONCLUSION 



In this contribution we study the dynamics of a con- 
densed atoms in an opto-mechanical system where it is 
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